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PRESSURE  DISTRIBUTION  ON  BLUNTED  FLAT  PLATES 
WITH  SURFACE  INCLINATION 


(Title  Unclassified) 


I  INTRODUCTION 

This  report  is  presented  as  a  partial  fulfillment  of  the  requirements 
of  Air  Force  Contract  No.  AF  33(6l6)-6692  dated  25  May  1959. 

The  object  of  this  report  is  to  present  and  analyze  the  experimental 
surface  pressures  measured  at  a  nominal  Mach  number  of  8  on  blunted 
flat  plates.  The  tests  are  conducted  in  the  B  tunnel  of  the  von  Karman  Gas 
Dynamics  Facility,  Arnold  Engineering  Development  Center.  Surface 
pressures  and  the  heat  transfer  to  the  model  were  measured  on  two  flat 
plates  blunted  by  a  semi-cylindrical  leading  edge.  The  nose  radii  of  the 
two  models  are  0.  10  and  0.  50  inches.  Both  models  had  a  chord  length  of 
1 1.  0  inches. 

In  this  report  the  experimental  pressure  distributions  are  presented. 
The  effects  of  angle  of  attack,  Reynolds  number  and  nose  radius  on  the 
surface  pressure  distribution  are  considered.  Tests  were  conducted  over 
a  range  of  0°  to  30°  surface  inclination.  The  free  stream  Reynolds  number 
varied  from  0.  43  x  I  0^  to  2.  00  x  10^  per  foot.  The  heat  transfer  data  will 


be  published  in  a  later  report. 
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Several  existing  theories  are  presented  and  their  comparisons  with 
experimental  data  are  made.  Also  presented  are  shadowgraphs  of  the 
flow  about  each  model  for  several  angles  of  attack. 
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empirical  (1/2)^^^  or  zero 

quantity  defined  by  Equation  43 

nose  drag  coefficient  defined  by  Equation  47 

nose  cylinder  diameter  -  inches 

stagnation  enthalpy  -  ft^/sec^ 

Mach  number 

pressure  -  Ibs/ft^ 

quantity  defined  by  Equation  41 

quantity  defined  by  Equation  46 

local  velocity  -  ft/sec 

nose  radius  -  inches 

Reynolds  number 

gas  constant  -  1716  ft^/sec^°R;  model  nose  radius,  inches 

entropy  -  ft^/  sec^;  arc  length  measured  from  stagnation  point,  inches 

temperature  °R 

velocities  defined  in  Figure  12 

chordwise  direction,  inches 

angle  of  attack  -  positive  nose  up 

ratio  of  specific  heats 


V 


GENERAL  APPLIED  SCIENCE  LABORATORIES,  INC. 


Page  4 

A  shock  stand-off  distance  (see  Figure  12) 

€  shock  wave  angle  measured  from  free  stream 

0  flow  deviation  measured  from  free  stream 

Mach  angle  (sin  *  -i-) 

M 

p  density  -  slugs/ft^ 

Subscripts 

GO  free  stream  conditions 

0  free  stream  stagnation  conditions 

Sg  stagnation  conditions  behind  a  normal  shock 

6  local  conditions  in  the  flow  field  outside  the  boundary  layer 

w  conditions  at  wall 
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A.  Tunnel 

All  tests  were  conducted  in  the  Mach  8.0,  50  inch  diameter  con¬ 
tinuous  axi symmetric  wind  tunnel  (Tunnel  B)  at  the  von  Karman  Gas 
Dynamics  Facility,  Arnold  Engineering  Development  Center,  Tullahoma, 
Tennessee.  A  complete  description  of  this  facility  may  be  found  in 
Reference  1. 

B.  Models 

Both  models  were  constructed  of  type  number  304  stainless  steel. 
The  surface  finish  was  nominally  50  microinches  while  the  skin  thickness 
was  0.25  inches.  The  nose  radii  of  the  two  models  are  0.10  inches  and 
0.  50  inches,  each  model  being  11.0  inches  long.  A  photograph  of  the  0.  50 
inch  radius  model,  before  being  positioned  in  the  tunnel,  is  shown  in 
Figure  1.  The  0,10  inch  radius  model  installed  in  the  tunnel  is  shown  in 
Figure  2.  The  model  sting  support  and  one  cooling  shoe  can  also  be  seen. 
Figure  3  shows  the  0.  50  inch  radius  model  with  the  side  and  top  plates 
removed  so  that  the  model  instrumentation  is  visible. 

Since  the  same  models  were  used  in  the  heat  transfer  tests,  slots 
were  milled  longitudinally  so  that  the  skin  thickness  in  the  vicinity  of  the 
thermocouples  was  nominally  0.050  inches.  The  thermocouples  were 
installed  in  these  slots.  A  more  thorough  and  comprehensive  description 
of  the  heat  transfer  aspect  of  the  models  will  be  given  in  subsequent  reports 
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Static  pressure  tubes  1/16"  O.D.  with  a  0.012"  wall  thickness 
were  installed  on  the  model.  These  orifices  were  placed  in  staggered 
rows  to  one  side  of  the  center  line.  Thermocouples  were  symmetrically 
placed  in  identical  locations  on  the  opposite  side  of  the  center  line.  The 
staggering  technique  was  used  to  obtain  a  small  axial  placement  of  tubes. 
The  location  and  number  of  the  tubes  used  for  each  model  is  shown  in 
Figures  4  and  5  and  Tables  1  and  2  for  the  0.  10  inch  radius  and  0.  50  inch 
radius  models,  respectively. 

C.  Test  Conditions 

All  tests  were  conducted  at  a  nominal  stagnation  temperature  of 
900°F.  Three  stagnation  pressure  levels,  namely,  100  psia,  240  psia  and 
480  psia,  were  chosen  to  determine  the  effect  of  Reynolds  number  variation 
on  the  pressure  and  heat  transfer  results.  Data  obtained  at  the  100  psi 
level  was  inconsistent  and  therefore  was  not  included  in  this  report.  The 
surface  inclination  of  the  plate  varied  from  0°  to  30°  in  2  l/2°  increments. 
However,  only  representative  angles  of  attack  of  0°,  5°,  10°,  15°,  20° 
and  30°  were  chosen  for  this  report. 

The  variations  in  test  conditions  influenced  the  tunnel  boundary 
layer  which  in  turn  affected  the  tunnel  Mach  number.  The  test  conditions 


are  summarized  below. 


GENERAL  APPLIED  SCIENCE  LABORATORIES.  INC. 


PAGt  7 


Model  11(A) 


Pjj  =  240  psi 

M  =  8.  03 

Pq  =  480  psi 

M  =  8.03 

Model  11(B) 

Pq  =  240  psi 

M  =  8.03 

Pq  =  480  psi 

M  =  8. 08 
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IV  DISCUSSION  OF  EXPERIMENTAL  RESULTS 

All  the  pressure  data  used  in  this  report  are  "preliminary  unchecked 
data".  Since  preliminary  data  were  used  certain  inconsistencies  in  the 
data  arose  and  were  unaccountable  for.  Therefore,  questionable  data  were 
usually  excluded  from  this  report. 

Because  of  the  small  nose  radius  on  Model  11(A),  the  nose  region 
(O^x/diO.S)  was  uninstrumented.  Results  of  the  pressure  distribution  on 
the  nose  region  of  Model  11(B)  are  shown  in  Figures  6  through  11.  Modified 
Newtonian  plus  Prandtl -Meyer  pressure  distribution  closely  approximates 
the  experimental  distribution. 

In  general,  the  data  are  self-consistent  and  comparison  with  several 


existing  theories  will  subsequently  be  discussed  in  this  report. 
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V  THEORETICAL  ANALYSES 

A.  Two-Dimensional  Rotational  Characteristics 

The  method  of  characteristics  for  supersonic  rotational  flow  can 
be  utilized  to  predict  the  flow  field  about  blunted  two-dimensional  bodies. 
However,  to  begin  the  computation,  properties  in  the  region  between  the 
shock  and  body  in  the  supersonic  region  must  be  known.  Thus,  the  computa¬ 
tion  procedure  that  was  utilized  was  divided  into  two  categories: 

1.  Subsonic- Transonic  Flow  About  a  Circular  Cylinder 

2.  2-D  Rotational  Characteristics  over  an  Arbitrary  Body  | 

Blunted  by  a  Circular  Cylinder 

A  description  of  each  category  with  the  pertinent  equations  follows. 

1.  Subsonic-Transonic  Region 

The  flow  field  about  two-dimensional  circular  cylinders  has 
been  analyzed  by  several  authors  (Refs.  2,  3  and  4).  In  each  method,  the 
technique  involved  in  computing  the  flow  field  is  quite  laborious  and, 
consequently,  time  consuming.  Since  the  object  of  this  analysis  is  to 
define  the  properties  in  the  transonic- super  sonic  region  of  the  flow  and  not 
the  details  in  the  subsonic  region,  an  approximate  scheme  was  devised. 

This  method  proceeds  by  assuming  a  shock  shape.  The  rela¬ 
tion  defining  the  shock  shape  must  be  "even”,  hence,  a  shock  of  the  following 
form  was  used  (Figure  12) 


=  'b  +  +  ‘■Z  *  ’>4 
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Since  the  shock  wave  radius  in  the  vicinity  of  the  stagnation  region  is  known 
to  be  the  best  approximated  by  a  concentric  cylinder  of  radius  (r^  +  A),  the 
square  term  must  be  eliminated.  Therefore,  a  shock  shape  of  the  form 

’^s  =  ^b  +  ^  +  ^4  (2) 

was  used,  where  the  parameter  b^  is  to  be  determined  from  mass  flow 

considerations.  That  is,  equating  the  mass  flow  entering  the  shock  (OA) 
to  the  mass  flow  passing  through  (AB)  (Figure  12)  provides  the  necessary 
cirterion  to  determine  b^. 

Equating  mass  flow,  and  retaining  only  the  linear  terms  in 
the  derivatives  one  obtains: 


i(A+b404) 


P  V  Ya  = 
^CD  00  A 


9p  9u 


i(A+b40^) 

Since  the  quantity  b^  can  be  expected  to  be  small,  only  linear  terms  in  b^ 
were  retained  throughout  the  entire  analysis.  The  pressure  at  point  B  was 
obtained  from  Newtonian  flow.  The  velocity  u,  and  density  p  were  then 
obtained  from  compressible  flow  relations. 

From  Crocco’s  theorem, 

-2i.  Grad  S  =  Curl  V  x  V  (4) 
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expanded  in  cylindrical  coordinates  one  obtains 
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iR  ar 


as  /  u  ,  9u  I  av 

—  =  -  u(- +  -^  -7  ae  ) 


But  at  the  body 


hence 


as  n  J  I  -  n 

dr  r  aO 


t'b 

Differentiating  the  relation 


P  =  P, 


E  -  'Tj'] 


in  the  r  direction,  one  obtains  (g~)^  expressed  in  terms  of  (  a^)g*  Hence 
the  first  integral  of  Equation  (3)  can  be  evaluated  in  terms  of  b^  > 


Defining 

,  T]  d  ro 
P  =  tan'*  i —  — — 
Ug  d0 


4b4  03 

tr  +  A+t'4  0'^ 


which  when  simplified  can  be  written  as 
P  =^1^4 


where 


"  rg+A 

In  order  to  evaluate  the  second  integral  of  Equation  (3),  the  velocity 


and  density  and  their  derivatives  immediately  behind  the  shock  must  be 
determined,  as  a  function  of  b4.  Within  the  approximation  considered,  the 
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velocity  components  behind  the  shock  can  now  be  expressed  in  the  following 
form: 


u  =  Uq  +  b4Uj 

V  =  Vo  +  b4  Vj 

P  =  Po  + 


(10) 

(11) 

(12) 


where 


u.  =  Sin  0 
o  cx> 


ui  - 


Pi 


(y  +  Uu 


CO 


cos  0  ^  2U*  SinJ.  - 

^  oo  0  cos 


-  2  U^j  Sin  0  Sin  20 


be[^ 


U,x)  +  Ujo  cos  2  0  -(V-I) 


I 


(y+1)Uoo 

Pi 


Sin  0  tan  0-|y  U,^  +  cos  20-  (^-1)  cos  0^(y-1)  + 2U^Sin2  0) 


r 

(> 

w 


V,  = 


‘  (y+1)  U 


GO 


Sin0|2U^  Sin^e-^^Q  (vU^+U^  cos  20  -  y  +  \^ 


+  2U^  cos  0  Sin  20 


n 


(Y  +  1)  cos^  0 

(Y-  1)  cos^  0  +  2 


oo 


2(y+l)M^  |3jSin20 
1^- 1)m4)  cos^  0  +  ^  ^  °° 

and 


U 


00 
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In  addition,  the  entropy  immediately  behind  the  shock  (Point  A, 
Figure  12)  can  b.e  expressed  as 


AS 


AS, 


AS, 


(■Y-I)  -  1)  +15^(7-  1) 


(13) 


where 


(\-l) 


AS, 


o  _ 


2yN^Cos^0-(y-l) 


log 

R  y+l 


-  \  log 


(y+l)Mj|j  Cos^  0 


and 


(y  -  1)  )  =  Pj  Sin  20 

R 


2yM2 


00 


(y-l)M^  Cos^  0 
y(y-l)l^  y 


2yM^Cos^0-(y-l)  (y -1)M^  Cos^  0+2  Cos^  0 


From  Equations  (10),  (ll)  and  (12)  the  derivatives  with  respect  to  0  of 
the  following  quantities;  u,  v  and ,  can  be  obtained  in  the  form 


d0  '  d0  ^  d0 


(14) 


With  the  system  of  equations 
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,,-4  )i^  mi-4  J  -  +  =  0 

r  9  0  v"*  or  r  r  av  or 


9S  _  .u  ^  ^  8v . 

VR  ar  dr  '  T  dO' 


a^  8S  _  .  u  8u 

9R  80  “  r 


r  80 


ds  8S  8S  dr 

dO  dO  87  ^ 


du  8u  8u  dr 
dO  ~  80  ^  8r  d  0 


dv  _  8  V  ^  dr 

dQ  8  0  8r  dO 


The  six  unknowns 

8u  9u  8v  8v  8^  8S 
W  ^  ’  ¥e  ’  9r  ’  90  ’  8r 

can  be  determined. 


The  remaining  term  to  be  evaluated,  ,  can  be  obtained  from 
the  continuity  relation 

13  3 

7  W  W  ^P'")  ~  ° 


That  is, 


p  8u  D  ^ 

rv  90  V  9r  rv  ^ 
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where  every  term  in  the  square  bracket  is  known  at  the  shocks  in  terms 
of  b4. 

The  second  integral  of  Equation  (3)  can  now  be  evaluated  in  terms 
of  b4.  Hence,  Equation  (3)  is  now  reduced  to  a  linear  algebraic  equation 
which  will  yield  the  value  of  b4. 

In  order  to  perform  the  numerical  computation,  it  is  necessary 
to  know  the  shock  detachment  distance  A.  The  results  of  the  constant 
density  solution  (References  2,  5  and  6)  yield  detachment  distances  which 
are  approximately  55%  too  small  at  Mach  3  and  45%  too  small  at  Mach  5. 
The  first  order  solution  of  Whitham  (Reference  2)  is  given  by 


Ti,  =  Ib  .  I  .  J_  log  ±L  (18) 

b  2\  ®  3 

and  X  =  Pz/P  1  '  ^ 

or  more  exactly  is  the  value  of  q  for  which  the  following  equation  is 
satisfied 


~  “!  I,(Xq)  r~  .  — 1  Ki(Xq) 

(X  +  1)K,(X)  -  \  K'i(X^  — —  -  ^  +1)  Ii(X)-X  Ii(xJ  ^ 

and  Kj  are  modified  Bessel  functions  of  the  first  and  second  kind, 
respectively.  Ijand  kJ  are  their  derivatives, 

ii^)  "  \  [^0^)  +  I2O0] 

K, 


(19) 


Ij  (X)=-  1  [Ko^)  +  K  (X) 


T1 

2 


(20) 


GENERAL  APPUED  SCIENCE  LABORATORIES.  INC 


pmi  1 6 


The  results  of  the  numerical  solution  of  Belotserkovskii  (Reference  3)  are 
compared  to  the  first  order  solution  of  Whitham  (Reference  2)  in  the  table 
below: 


Belotserkovskii 

Whitham 

(Ref.  3) 

(Ref.  2) 

M  «=  3 

A/rb  =0.703 

A/rb  " 

M  =  4 

=  0.  546 

=  0. 279 

II 

Ln 

=  0.481 

=  0.  265 

In  order  to  find  the  shock  detachment  distance  at  Mach  8.00,  the  first  order 
solution  and  the  "exact"  solution  of  Reference  2  (Eqs,  18  and  19)  were  solved 
and  the  results  were  A/ =  0.234  and  0.  279,  respectively.  These  results 
were  obviously  too  small.  To  determine  the  detachment  distance,  sim:e 
data  at  M  =  8.  0  was  not  available  in  literature.  Equation  (18)  was  used  as 
a  curve  fit  for  Belotserkovskii' s  data.  That  is 


where 


^  I  .  — 1 - ^  log  ^(\o  +  A\) 

2(ko+AK)  3 


(-L)  -I  =  4.000 

M  =  5.00 


AX=(\)  -(X)  =0.565 

M  =  8.00  M  =  5.00 
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The  value  for  the  detachment  distance,  A,  determined  from  the 
above  relation  is 

A/rjj  =  0.422 

The  corresponding  value  of  obtained  from  Equation  (3)  is 
b4  =  0.1815 

Thus  the  shock  shape  is  defined  as 

I®.  =  1.422  +  0.  1815  0^  (22) 

^b 

and  at  0  =  55° 

—  =  1.576 
^b 

With  the  shock  specified,  the  velocities  u  and  v  and  the  entropy 
S/R  and  their  derivatives  in  the  r  direction  are  known  at  the  body  and  at 
the  shock.  To  obtain  data  in  the  region  between  the  shock  and  the  body,  a 
third  order  equation  was  used.  In  each  case,  the  four  boundary  conditions 
were  used  to  solve  for  the  four  coefficients.  Hence  along  the  line,  0  =  56° 
the  following  properties  were  determined  at  several  points: 

x/rjj,  y/r]3,  S/R,  P/Pco’  ^ 

These  data  are  the  input  to  the  characteristics  analysis. 

2.  Two-Dimensional  Rotational  Characteristics 

The  analysis  of  the  flow  over  blunted  two-dimensional  bodies  has 
been  computed  by  the  method  of  characteristics  for  a  rotational  flow. 
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The  equation  of  motion,  when  the  total  enthalpy  Hq  is  constant 
in  the  flow,  can  be  expressed  along  the  characteristic  lines  in  the  following 
form 


1 

V  P" 


+ 


d8 

sin  p  cos  p 


=  0 


(23) 


along  the  characteristic  line 

tan  (6  +  p)  =  ^  (24) 


and  1  ^  - 

de 

(25) 

V  P 

sin  p  cos  p 

along  the  line 

tan  (0  - 

dx 

(26) 

and 

dS  =  0 

(27) 

along  the  streamline  defined  by 

=  tan  0  (28) 

dx 

If  the  velocity  vector,  total  enthalpy,  and  entropy  are  known  at  two  points 
A  and  B  (see  sketch),  the  position  of  a  point  C  at  the  intersection  of  the 
tangents  AC  to  the  characteristic  line  of  the  first  family  from  A  and  the 
tangent  BC  to  the  characteristic  line  of  the  second  family  from  B  can  be 
obtained  since  the  value  of  p  can  be  calculated  at  A  and  B. 

Assuming  that  the  tangents  AC  and  BC  can  be  substituted  for  the 


characteristic  lines  between  AC  and  BC,  the  first  approximation  to  the 
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variation  of  flow  properties  between  A 


and  C  and  B  and  C  can  be  obtained  from 


Equations  (23)  and  (25)  by  assuming  that 
the  coefficients  are  constant  and  equal 
to  the  corresponding  values  at  A  and  B. 


However,  the  variation  of  the  entropy 


along  the  characteristic  lines  must  be  known  in  order  to  determine  P  and  0 
at  point  C. 

The  assumption  that  the  entropy  is  a  continuous  function  in  the 

» 

region  of  the  flow  between  the  streamlines  that  pass  through  A  and  B,  plus 
the  fact  that  the  entropy  remains  constant  along  a  streamline  leads  to  the 


equations 


(AS),„  =  Sc-S, 


(^®)bc  ~  ■  ^b 


Sc  =  Sa  ' ac"^"’ac  “  ® a  "  <  dJ  )ac  '*0 


_  — n 

1^08  (0  + 


Sc  =Sb  t  (^)bc''^”>bc  CoI(C-  I.)  g 


where  ^  is  the  gradient  normal  to  the  streamlines,  and  An  is  the  normal 


distance  between  the  streamline  at  A  or  B  and  the  streamline  passing  at  C. 
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Therefore  assuming  that 

( ) 

Hn  ^  BC 


(M) 

dn  AB 


Sq  can  be  determined  once 


dS 


is  obtained  from  the  expression 


(32) 


-  Sa 


AB 


(33) 


X-  - 


C  -  ^A 


|Cos(0  +  pQ 


Sin 


[Cos{0  - 


The  pressure  P^;  and  the  flow  deviation  Oq  can  then  be  obtained  from  the 
relations 


Pc  = 


Qa-Qb^^  |iin  iiaCosI^A  I^B 

>^B  * 

V  Pb 


Si  n  Cos  M’  A  ^  B 


®  =  ®A  +  V  ^"A 


Pa  -  Pc 


Pa 


(34) 


(35) 


Shock  Conditions 

The  determination  of  points  on  the  shock  were  obtained  from  the 
following  shock  relations. 


(V  -  l)Mjjj  Sin^  €  +  2 


(y  +  1)  Mj  Sin  ^  e 

Sin  e 

(36) 

=  Cos  « 

(37) 

+  ^Tc 

(38) 

8 
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P*ce  2  I 


oo 


2  yMcx)  ‘  -  (\-l) 

(Y  +  1) 


and  the  characteristic  relation 


(39) 


1  Pa  -  Pc  ^  Qa-  Oc 

V  Sin^^Cosfi^ 


(40) 


The  conservation  of  tangential  momentum  across  the  shock  can  be  expressed 
as 

Vqo  =  V  Cos  0  +  tan  €  Sin  0  (41) 

The  properties  at  point  C  (see  sketch)  on  the  shock  must  satisfy  both  the 

shock  and  characteristic  equations.  Since  the  problem  was  to  be  performed 

on  an  automatic  computing  machine,  an  iterative  scheme  on  Equation  (41) 

was  decided  upon  as  being  the  most  expeditious  way  of  obtaining  a  solution. 

(1)  (2  ) 

The  shock  inclination,  ,  at  C  was  selected  and  then  was  determined 

from  Equation  (41). 


I 
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Boundary  Conditjons 


To  determine  properties  on  the  body 
(e.g.,  Point  C),  the  second  family 
characteristic  equation  (Eq.  25)  and  the 
and  the  geometry  of  the  body  are  used. 


Since  the  body  is  a  streamline,  Sa  “ 
The  flow  deviation  is  known  for  the  geometry,  thus  one  can  solve 


Equation  25  directly  for  P.  That  is, 


Pc  =  P] 


1  +  V 


®B  ■  ®C 


Sin  pb  cos  pg 


(42) 


The  above  equations  were  programmed  for  the  Bendix  G-15 
digital  computer.  In  order  to  systematize  the  computations,  the  flow  field 
was  divided  into  three  regions  (Figure  13). 

1.  Region  A  is  defined  by  the  second  family  line  from  point  S  and 
the  first  family  line  from  point  B. 

2.  Region  B  is  defined  by  the  body  and  the  second  family  line 
from  point  C. 

3.  Region  C  is  the  remaining  flow  field. 

The  computation  in  region  C  was  performed  along  second  family 
lines.  In  this  manner,  the  program  is  able  to  sense  terminal  points  of  the 


body  profile. 
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Attention  is  currently  being  devoted  to  minimizing  the  computer 
production  time  while  retaining  reliability.  The  problem  is  hampered 
solely  due  to  the  limited  capacity  of  the  computer  availabile.  However, 
modifications  are  .  eing  made  which  will  substantially  reduce  the  computer 
production  time.  Presented  in  this  report  are  the  characteristic  analyses 
for  two  angles  of  attack,  namely,  a  ==  0°  and  a  =  5°.  In  a  subsequent  report 
several  additional  cases  will  be  presented,  including  a  displacement  thick¬ 
ness  analysis  for  which  the  inviscid  two-dimensional  characteristic  will  be 
computed. 

B.  Modified  Newtonian  Pressure  Distribution 

This  method  for  predicting  the  pressure  distribution  on  the  surface 
of  a  blunted  flat  plate  is  discussed  in  detail  in  References  7  and  8.  On  the 
flat  portion  of  a  blunted  plate,  the  normal  to  the  free  stream  direction  is 
equal  to  the  compliment  of  the  angle  of  incidence  of  the  plate. 

The  equation  relating  the  pressure  to  the  plate  incidence  (Ref.  8) 
can  be  written  as 

P„/P„  =  Sin^a  +■  ^  Cos^  a  (43) 

e  p 

®e 

For  a  given  free  stream  Mach  number,  the  pressure  distribution  is  solely 
a  function  of  surface  inclination.  Therefore  on  the  constant  inclination 
surface  of  the  flat  plate,  the  pressure  as  predicted  by  Newtonian  theory  is 


constant. 
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C.  Modified  Newtonian  Plus  Prandtl -Meyer  Pressure  Distribution 

This  method  attempts  to  modify  the  body  pressure  distribution 

predicted  by  the  Newtonian  theory  from  the  point  on  the  body  surface  where 
the  pressure  and  pressure  gradient  is  equal  to  that  predicted  by  Prandtl- 
Meyer  flow.  A  more  complete  description  of  its  use  is  given  in  Ref.  8. 

The  Prandtl -Meyer  correction  to  Newtonian  theory  tends  to 
alleviate  somewhat  the  overexpansion  that  the  modified  Newtonian  theory 
predicts. 

D.  Wedge  Solution 

This  method  is  simply  the  solution  to  a  pointed  two-dimensional 
wedge  inclined  to  the  free  stream  at  the  same  angle  as  the  blunted  flat 
plate.  Theoretically,  at  a  large  distance  downstream  from  the  nose,  the 
pressures  predicted  by  the  blunted  flat  plate  and  the  wedge  solution  should 
be  identical.  For  a  given  length  flat  plate,  the  wedge  technique  would 
obviously  be  more  applicable  for  a  smaller  nose  radius  on  the  blunted  plate. 

E.  Creager*s  Approximate  Pressure  Distribution  for  Blunted 

Flat  Plates 

This  technique,  as  outlined  more  completely  in  References  9  and 
10,  assumes  that; 

1.  A  detached  bow  shock  wave  exists 


2.  The  flow  field  is  comprised  of  a  viscous  boundary  layer 
adjacent  to  the  surface  of  the  body,  plus  an  inviscid  region 
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between  the  boundary  layer  and  bow  shock. 


p*«(  2  5 


3.  The  displacement  effect  of  the  boundary  layer  perturbs  the 
inviscid  flow  field. 

If  the  above  assumptions  are  valid  then  the  perturbation  may  be 


expressed  as: 


where 


0.  865  '^w 

^6=  lip"  17  +  (V-  1)  y 


From  References  2,  3,  4,  and  7  of  Reference  9  it  is  further  assumed  that 
the  inviscid  pressure  is  given  by  a  sum  of  two  terms:  one  including  a 
contribution  due  to  the  surface  inclination  and  the  second  due  to  the  shock 
wave  shape. 

Therefore  it  follows  that 

p  _  ^  P 

p  (46) 

^oo  ^00  oo 


Combining  Equations  (42)  and  (44)  and  multiplying  both  sides  by  Pqq/Ps 


Pse  Lfse  L 


The  inviscid  pressure  terms  (Eq.  (46))  merit  some  further  consideration: 


is  assumed  to  be  a  function  of  the  surface  inclination. 
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As  mentioned  in  Reference  9,  this  term  is  assumed  to  be  given 
by  the  modified  Newtonian  impact  theory.  However,  as  pointed  out  in 
Section  VD  of  this  report,  the  surface  pressure  should  approach  that  value 
given  by  an  equivalent  wedge. 

It  is  for  these  reasons  that  the  chordwise  pressure  distributions 
as  given  by  Reference  9,  have  been  calculated  by  assuming  a  modified 
Newtonian  distribution  for  the  large  nose  radius  model  and  an  equivalent 
wedge  value  for  the  small  nose  radius  model. 

The  second  contribution  to  the  inviscid  surface  pressure  (that  due 
to  shock  wave  shape)  has  been  computed  by  utilizing  the  blast-wave  theory. 
From  Reference  9  this  is  assumed  to  be  given  by; 

(48) 

where 


Pf: 


^  BC 


oo 


'D 


X  +  A 


=  0.1  12,  and  Cj-j  is  given  by: 


^  P  M  '  ^ 
2  oo  ^oo  2 


P  tan  0  dx 
w 


(49) 


For  the  body  at  angle  of  attack,  the  inclination  of  the  circular  nose  is: 

0  =  90  -  a 

w 

Also 


X 


r 


r  cos  Q 
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The  limits  of  integration  for  Equation  (47)  have  been  modified  in  order  to 

incorporate  the  portion  of  the  nose  which  most  influences  the  drag  coefficient. 

Therefore  Cj)  as  suggested  in  Reference  9  was  modified  to: 

r(l- sin  a) 


^oo  ^oo 


v; 


P  tan  0^  dx 


(50) 


r(l-cos  q) 


The  contribution  to  the  pressure  for  the  zero  angle  of  attack  case  of  each 
term  in  Equation  (45)  can  be  seen  in  Figures  14  and  15  for  each  of  the  two 
models. 


GENERAL  APPUED  SCIENCE  LABORATORIES,  INC 


Pa«  28 

VI  DISCUSSION  OF  RESULTS 

The  pressure  distribution  on  two  blunted  flat  plates  were  experimental¬ 
ly  determined  at  a  nominal  Mach  number  of  8.  0  for  a  0°  to  30°  angle  of 
attack  range.  Tests  were  conducted  over  a  free  streeum  Reynolds  number 
range  of  0,43  x  10^  to  2.00  x  10^  per  foot.  The  two  models  were  blunted 
by  semicylindrical  leading  edges  with  nose  radii  of  0.  10  inches  and  0.  50 
inches.  Both  models  had  a  chord  length  of  I  1.0  inches. 

The  experimental  pressure  distribution  on  the  cylindrical  nose  of 
Model  11(B)  (l/2"r)  has  been  plotted  and  compared  to  the  modified 
Newtonian  impact  theory.  The  results  are  shown  for  several  angles  of 
attack,  namely,  0°,  5°,  10°,  15°,  20°  and  30°.  It  is  generally  concluded 
that  the  modified  Newtonian  theory  is  in  good  agreement  with  experiment 
in  the  nose  region  of  blunt  bodies  (Ref.  6,  page  257).  One  will  note  from 
the  figures  (6  to  1  1)  that  the  agreement  with  experiment  is  not  as  good  as 
can  be  expected.  It  appears  that  the  model  angle  of  attack  is  greater  than 
the  value  indicated  on  the  "Preliminary  Data  Sheets".  For  example,  if  the 
model  angle  of  attack  was  18°  or  19°  instead  of  15°  (Fig.  9),  the  modified 
Newtonian  theory  would  be  in  good  agreement  with  experiment. 

The  pressure  distribution  on  the  1/2"  nose  radius  flat  plate  model 
tends  to  approach  the  Newtonian  value  for  all  angles  of  attack  (Figs.  1 6  to  21). 
In  addition,  for  the  low  angle  of  attack  range,  blast  theory  plus  the  viscous 
correction  (Refs.  9  and  10)  accurately  predicts  the  pressure  variation  on 
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the  flat  plate.  The  agreement  between  the  analysis  of  Reference  9  in  the 
forward  region  of  the  flat  plate  is  quite  poor  for  larger  angles  of  attack. 

The  asymptotic  solution  for  Creager's  analysis  for  the  l/2"r  model  is  the 
Newtonian  value.  Therefore,  even  for  the  larger  angles  of  attack,  the  blast 
theory  predicts  reasonably  well  the  pressures  on  the  rearward  portion  of 
th^  plate.  The  rotational  characteristics  analysis  has  been  computed  for 
a  a  0°'and  5°  and  the  results  agree  well  with  experiment. 

j 

iFcir  the  0.  10"  nose  radius  model,  the  pressure  distribution  on  the  flat 
plate  tends  to  approach  the  wedge  value.  Hence,  the  wedge  value  was  used 
as  the  asymptotic  solution  in  the  analysis  of  Reference  9.  The  blast  solution 
plus  the  viscous  correction  accurately  predicts  the  pressure  distribution  on 
the  flat  plate  over  the  entire  angle  of  attack  range  considered  (Figs.  22  to  27) 
For  comparison,  the  zero  blast  solution  (B  *  0,  Eq.  (48))  is  also  shown. 

For  higher  angles  of  attack  (20°-30°)  this  solution  appears  to  be  in  slightly 
better  agreerrient  than  the  blast  theory.  This  is  to  be  expected,  since  the 
majority  of  streamlines  that  "wet"  the  body  pass  through  the  oblique  portion 
of  the  bow  shock.  Hence  the  blunted  flat  plate  acts  more  like  a  sharp  wedge. 
Several  shadowgraphs  of  each  model  are  shown  in  Figures  28  to  39  for  the 
entire  angle  of  attack  spectrum. 
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Figure  3  -  0.  50  Inch  Nose  Radius  Flat  Plate  Model  Instrumentation 
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Figure  7  -  Pressure  Distribution  on  the  Nose  Cylinder  of  the  0.  50  Inch  Nose  Radius 
Flat  Plate  Model  -  a  =  5° 


Figure  9  -  Pressure  Distribution  on  the  Nose  Cylinder  of  the  0.  50  Inch  Nose 
Radius  Flat  Plate  Model  -  a  =  15° 


Figure  10  -  Pressure  Distribution  on  the  Nose  Cylinder  of  the  0.50  Inch  Nose 
Radius  Flat  Plate  Model  -  a  =  Z0° 


Figure  16  -  Axial  Pressure  Distribution  for  the  0.  50  Inch  Nose  Radius  Flate  Plate  Model 


Figure  18  -  Axial  Pressure  Distribution  for  the  0,50  Inch  Nose  Radius  Flat  Plate  Model 


Figure  19  -  Axial  Pressure  Distribution  for  the  0.50  Inch  Nose  Radius  Flat  Plate  Model 
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Figure  20  -  Axial  Pressure  Distribution  for  the  0.  50  Inch  Nose  Radius  Flat  Plate  Model 
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Figure  26  -  Axial  Pressure  Distribution  on  the  0.  10  Inch  Nose  Radius  Flat  Plate  Model  -  a  =  20 
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Figure  28  -  ShadJ^graph  of  the  Flow  Over  the  0.  50  Inch  Nose  Radius 
Flatmate  Model 
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